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Abstract 

We study the contribution of new sets of two-loop Barr-Zee type diagrams to the anomalous 
magnetic moment of the muon within the two-Higgs-doublet model framework. We show that 
some of these contributions can be quite sizeable for a large region of the parameter space and 
can significantly reduce, and in some cases even explain, the discrepancy between the theoretical 
prediction and the experimentally measured value of this observable. Analytical expressions are 
given for all the calculations performed in this work. 
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1 Introduction 


Now that a SM-like Higgs particle has been experimentally discovered Si, the possibility of an 
enlarged scalar sector becomes very plausible. In this analysis we are going to use the anomalous 
magnetic moment of the muon as a probe for new physics and study new contributions to this ob¬ 
servable within the two-Higgs-doublet model (2HDM) framework. The anomalous magnetic moment 
of the muon has been extensively analysed within the Standard Model (SM) and its numerous exten¬ 
sions. Even if the SM prediction still suffers from large theoretical uncertainties (mostly hadronic and 
electroweak) it is a nice place to look for new physics. The latest result for the discrepancy between 
the SM prediction and the experimental measured value is given by [6 

A= a e * p - a® M = 262(85) x HT 11 . (1) 


28 


Here we will study the two-loop Barr-Zee type |29| contributions to Athat have not been anal¬ 
ysed previously within the 2HDM. We show that some of these diagrams can bring rather sizeable 
contributions for a quite large region of the parameter space and therefore can reduce the value of 
the difference between theory and experiment given by ©■ We also show that other sets of these 
type of diagrams bring small contributions and can be safely discarded. For the calculations we use 
the most generic Higgs potential and the generic Yukawa structure of the aligned two-Higgs-doublet 
model (A2HDM) 1 30 1. Thus, we also re-examine the classical Barr-Zee type diagrams |6][7J[31-42 
expressing their contributions in terms of the three independent complex alignment parameters s u ,d,l- 
All the results are given in analytical form. The phenomenological analysis is made assuming a CP- 
conserving Lagrangian. However, all the generic formulae given in this work can be used for future, 
and more complete, analyses without assuming CP-conservation. Additional constraints coming from 


the flavour sector and global fits to the LHC data are also taken into account 43-52 


In the first part of this paper, section[2j we present the relevant features of the A2HDM. In section[3] 
we present the one-loop results in terms of the generic A2HDM parameters. In section [4] we present 
the classical two-loop Barr-Zee results and the calculation of the new sets of this type of diagrams 
that can potentially bring sizeable contributions to the anomalous magnetic moment. Section [5] is 
dedicated to the phenomenological analysis for the CP-conserving case and the presentation of the 
relevant contributions. Finally, we conclude in section [6] with a brief summary of our results. One 
appendix is also given, with technical details for the calculation of a particular set of Barr-Zee type 
diagrams. 


2 







$1 = 


$0 = 


(2) 


2 The Aligned Two-Higgs-Doublet Model 

The 2HDM extends the SM with a second scalar doublet of hypercharge Y = ^. It is convenient 
to work in the so-called Higgs basis (<3?i, $ 2 ), where only one doublet acquires a vacuum expectation 
value: 

G + 

^ (v + 5i + iG°) 

where Cr and G° denote the Goldstone fields. Thus, $1 plays the role of the SM scalar doublet with 
v = (y/2 Gf)^ 1 / 2 = 246 GeV. The physical scalar spectrum contains five degrees of freedom: two 
charged fields H ± {x) and three neutral scalars x ) = {h(x), H(x), A(x)}, which are related with the 
Si fields through an orthogonal transformation <Pi(x) = IZijSj(x). The form of the 1Z matrix is fixed 
by the scalar potential, which determines the neutral scalar mass matrix and the corresponding mass 


H+ 

[^(S 2 + iS 3 ) J 


eigenstates. A detailed discussion is given in 143 -45 . In general, the CP-odd component S 3 mixes 


with the CP-even fields 5^2 and the resulting mass eigenstates do not have a definite CP quantum 
number. If the scalar potential is CP symmetric this admixture disappears; in this particular case, 
A(x) = S 3 (x) and 




( 3 ) 


cos a sin a 
— sin a cos a 

Performing a phase redefinition of the neutral CP-even fields, we can fix the sign of sind. In this work 
we adopt the conventions < M R and 0 < a < n, so that sin a is positive. 

The most generic Yukawa Lagrangian with the SM fermionic content gives rise to FCNCs because 
the fermionic couplings of the two scalar doublets cannot be simultaneously diagonalized in flavour 
space. The non-diagonal neutral couplings can be eliminated by requiring the alignment in flavour 


space of the Yukawa matrices 30 ; i.e., the two Yukawa matrices coupling to a given type of right¬ 


handed fermions are assumed to be proportional to each other and can, therefore, be diagonalized 
simultaneously. The three proportionality parameters (/ = u, d, l ) are arbitrary complex numbers 
and introduce new sources of CP violation. 

In terms of the fermion mass-eigenstate fields, the Yukawa interactions of the A2HDM read (30| 

V2 


C Y = ~H+{u VM d V R - MlVV L ] d + Q vM{P R l\ 


1 


£ vfrf [.fM f V R f\ + h.c. 


( 4 ) 


vU 


where V Ri l = lj ^ 75 are the right-handed and left-handed chirality projectors, Mf the diagonal fermion 
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mass matrices and the couplings of the neutral scalar fields are given by: 


Vd'l — Tin + (^i2 + i ^ 3 ) Sd,i > 


2/u* = Un + ( Ui2 - i 7 ^ 3 ) 


( 5 ) 


The usual models with natural flavour conservation, based on discrete Z 2 symmetries, are recovered 


for particular (real) values of the couplings 30 . The coupling of a single neutral scalar with a pair 
of gauge bosons takes the form (V' = W, Z) 


„ _ -t) „SM 

9tp°VV — / '-il 9hVV ! 


( 6 ) 


which implies gfyy + gjjyy + 9avv = Owv) 2 - Thus, the strength of the SM Higgs interaction i 


is 


shared by the three 2HDM neutral bosons. In the CP-conserving limit, the CP-odd field decouples 
while the strength of the h and H interactions is governed by the corresponding cos a and sin a factors. 
Again, for further details about the interaction Lagrangian as well as the Higgs potential, needed for 


the calculations in this work, see 143 - 45 . 


3 One-loop contribution 


At the one-loop level, the contribution of the 2HDM extension of the SM to the anomalous magnetic 
moment of the muon is given by the two well-known diagrams shown in Fig. [lj The explicit expressions 
for these contributions, in terms of the most generic Higgs potential and the A2HDM Yukawa structure, 
are given by 


, , m.. 

A a (“) = _IF 

" 8 tt 2 u 2 


E 


m, 


M\ 


Vi L 


n , v °x 2 f 1 , x 2 (2 — X ) 

81V ‘ L 1 (mJ/Mjo)i 2 - x + 1 


+ ImW? ) 2 7„ M 


^ —x 3 


for the neutral Higgses and 


m. 


Aa (b) = _ 

^ 8 tt 2 u 2 Ml 


H± 


I 12 f 1 1 , ® 2 (1 - x) 

lQl Jo (ml/M% ± )x(l-x)-x' 


(7) 


( 8 ) 


for the charged Higgs. These contributions have been previously analysed in 6,31, 36,[53]-56|. 

It’s a known fact that the two-loop Bar-Zee type diagrams dominate over the one-loop contribu¬ 
tions. The two loop contributions have a loop suppression factor of {ol/tt) but also have an enhancement 
factor of (M 2 /m 2 ), where M stands for the mass of heavy particles running in one of the loops: M H ±, 
rrit■ 0 , etc. This last factor usually dominates over the first one. Furthermore, in the usual Z 2 
models, there is an extra enhancement (suppression) factor from tan [3 (cot (3) for some diagrams. In 
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Figure 1: One-loop contribution to A in two-Higgs-doublet models. 

the aligned model there is a lot more freedom to independently enhance or suppress any contribution 
through the alignment parameters We shall see next, that for somewhat large values of these 
parameters, there are new Barr-Zee contributions that have never been taken into account, and can 
bring quite sizeable contributions to (g — 2)^. 



4 Two-loop contribution 


The Barr-Zee type contributions with an internal photon, he., Fig. [2] diagrams (1) and (2), have 
been extensively analysed within the 2HDM and also in minimal super-symmetry (MSSM) framework 
31-42 . Diagram (3) from Fig. [2] is also of the Barr-Zee type and could, in principle bring 


important contributions. Given that the coupling to a pair of gauge bosons of the recently discovered 
scalar particle is close to the SM prediction 143], one expects the contributions from the remaining 
scalars to be somewhat suppressed (by a factor 1Z, \). However, we shall see that this statement is not 
correct, and that this contribution is quite sizeable. 

Similar contributions to the ones shown in Fig. [2] but with the internal photon replaced by a Z 
boson have been also analysed in the literature [33]. These contributions have a relative suppression 
factor of order 1CP 1 2 . This factor is in part due to the vectorial couplings of Z to leptons, which 
are the only ones that survive for both scalar and pseudo-scalar bosons [33], and in part from the Z 
propagator which introduces a new mass scale Mz- Therefore we will ignore these contributions in 
our present analysis. 

This is, pretty much, the summary of all the mechanisms that are usually considered in the 
literature. However, there is no reason a priori to discard other similar Barr-Zee contributions with 
a charged Higgs H± substituting the neutral scalars tp®, and a W boson substituting the internal 
photorQ These diagrams are illustrated in Fig. [ 3 ] On one hand, one expects a relative suppression 


1 Similar contributions, however, with sfermionic loops within the MSSM framework have been previously analysed 


in 57 . 
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Figure 2: 


Two-loop Barr-Zee type (with an internal photon) contribution to A in two-Higgs-doublet models . 





Figure 3: Two-loop Barr-Zee type (with a charged Higgs and an internal W boson) contribution to A in 
two-Higgs-doublet models. 


factor with respect to the contributions of the diagrams from Fig. [2] due to the propagator of the 
W boson (note that in this case we don’t have the additional suppression factor due to the gauge 
boson couplings to leptons, as in the Z case). On the other hand, one must also expect to be able 
to re-enhance these contributions with the (or tan f3) parameters, and therefore, obtain sizeable 
contributions at least in some regions of the parameter space. 

In this analysis we shall calculate the contribution from these new diagrams and demonstrate, 
that in fact, all of these new sets can bring rather sizeable contributions to the anomalous magnetic 
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Figure 4: Generic two-loop Barr-Zee type contributions, with two internal charged Higges (left) and two internal 
W bosons (right), to A in two-Higgs-doublet models. 

moment of the muon in a quite large region of the parameter space. For completeness we shall also 
present the classical two-loop results in terms of the most generic Higgs potential and in terms of the 
generic Yukawa texture of the A2HDM. 

Before moving on to the next section and presenting the analysis, there are a couple of related 
cases that are worth discussing. They are shown in Fig. [4j where the grey circles stand for the same 
loop contributions as in Fig. [3] (excluding the fermionic loops for diagram (B) which is just a pure 
SM contribution). The contribution from the first case (A), will have a relative suppression factor 
rrijjM^y with respect to the contributions of diagrams from Fig. |3] so we can safely discard it. The 
contribution coming from the second set, Fig. [4] (B), does not have this suppression factor, thus we 
can expect, at least in principle, a rather sizeable effect. Details of the the full calculation of this last 
set of diagrams, together with other technical details are given in appendix |A| Roughly one obtains a 
contribution of 0(10 X1 ) which is rather small and we shall not include it in this analysis. 

Next we move on to the analysis of the set of diagrams shown in Fig. [3] which is the main goal of 
our paper. 

4.1 Gauge invariant effective vertices 

The calculation of the two-loop Barr-Zee type diagrams can be separated in two parts. We will first 
calculate the ip® — 77 and H + — 7 W + one-loop effective vertices and obtain analytical and rather 
simple expressions. With these expressions, the calculation of the second loop becomes quite trivial. 
The effective vertices can be written in a generic gauge-invariant transverse form: 

iT^ = i (g^k • q - Wq u ) S + i k a q p S , (9) 
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Figure 5: Feynman rule for the gauge-invariant one loop effective vertices — 77 and H + — W + 7 . 

where q 11 ' is the momentum of the incoming real photon and k v is the momentum of the out-going 
virtual gauge boson (see Fig. [HJ, and where S and S are scalar form factors. In order to obtain 
this expression we have considered the most generic Lorentz structure for the vertex, and we 
have imposed the electromagnetic current conservation qy T^ v = 0. All terms proportional to q^ have 
also been eliminated as they cancel when contracted with the polarization vector of the photon. As 
the W boson is off-shell, in the actual calculation of the effective vertex there will also appear some 
other Lorentz structures than the ones shown in ©• However in some cases, these gauge-dependent 
contributions vanish when calculating the second loop or they are cancelled by some other non Barr- 


Zee terms, as it is nicely shown in 58 . If this was not the case, when summing the proper non Barr-Zee 


contributions to the gauge dependent Barr-Zee terms, the result must be gauge independent. As the 
gauge dependence from the Barr-Zee terms is cancelled by other sub-dominant topologies, we also 
expect this contribution to be sub-dominant. Therefore, we shall discard these terms in our analysis. 

The gauge independent contribution from each set represented by the generic topologies in Fig. [2] 
and Fig. [ 3 ] is transverse by itself, i.e., of the form given in ([ 9 ]); we can therefore decompose the results 
into eight separate contributions. For the tp { - — 77 effective vertex S = S^y + S^+S^ and S = Syy, as 
for the H + — yW + vertex we have S = S^ + S^ + and S = S^y Note that the only contributions 
to the k a qp structure come from the fermionic loops. Furthermore, one can adopt our strategy 
from 45 , and further simplify the calculations of Sy) by only considering the terms that contribute 


to the structure k fl q u . 

It is worth mentioning the following technical detail. When performing the calculations for the first 
loop, after introducing the Feynman parametrization and after integrating over the four-momentum, 
one obtains a denominator similar to 


[k 2 x(x — 1) + M 2 x + M 2 ( 1 — x) + k ■ q2yx(l — x)] 


-1 


( 10 ) 


where M a ^ are the masses of heavy particles running in the loop, i.e., Mw, mt, M H ± , etc. It is 
a very common assumption that the photon is “soft” so one can ignore the k ■ q term as a good 
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approximation. This term, in fact, can be safely ignored without making any assumptions on the 
“softness” of the photon. Keeping track of this term, one can observe that it simply vanishes when 
calculating the second loop integral. However, this happens accidentally for diagrams (1) to ( 6 ); for 
the WW'y effective vertices calculated in appendix [A] this is not always the case. Thus, having checked 
that these terms play no role in our present case, we will discard them already at the one-loop level in 
order to give simpler and more elegant expressions for the form factors SA and . After performing 
the four-momentum loop integral we obtain the following expressions for the scalar form-functions 


a m 


S = y ^Q2 N f Re ^, f d 2X(1-X)-1 

y ) nv f c v / ; k 2 x( 1 — x) — m 2 f 


(ii) 


a m 


S(d = -- Q 2 tNL Im(iy ! f' 1 ) / dx ---~ 

7 tv * c K J Jo k 2 x(l — x) — m 2 f 


hf 


5 (2) “ Z TT X V°H+H- 


dx 


x{x — 1 ) 


2 vr ' n n J o k 2 x( 1 - x) - M 2 h± ’ 

for the tp® — 77 vertices with a fermionic or a charged Higgs loop, in agreement with 
third diagram, we find 


41 


( 12 ) 

(13) 
As for the 


73) = £ 


7 T— Kil I dx 
2 nv Jo 


i M 2 v x(3x(4x — 1) + 10) — M 2 o x(l — x) 




k 2 x(l — x) — M 2 V 

The new gauge-invariant scalar form factors coming from diagrams (4) to ( 6 ) are given by: 


(14) 


%) = 


aN c \V tb \ 2 r 1 ^ \Q t x + Q b ( 1 - x)\ [s u m$x 2 - <^mg(l - x) 2 ] 


27T V s v 


k 2 x( 1 — x) — m 2 ( 1 — x) — m 2 x 


- . aN c \V tb \ 2 [ l [Qtx + Q b (l-x)\[-^ u rnfx + ^ d ml(x-l)\ 

S ( 4 ) = l --- / dx - , 0 -r- Try- -r-- 

2ttv s w Jo k z x{l — x) — m£(l — x) — m^x 


<%) = 


a 


47r v s v 


Jo 


{M 2 h± + Myy - M 2 ){ 1 -x)- 4 M 2 , 

• i 


dx T 

o ' k 2 x(l - x) - M^x - M 2 0 (l - x) 


Vi 




a v 
47 


(^-*2 — *^ 3 ) 


dx 


x 2 (x — 1 ) 


0 A; 2 x(l — x) — M 2 ± x — M 2 0 (1 — x) ’ 


(15) 

(16) 

(17) 

(18) 


with s w = sin 0 W , and 0 W the weak mixing angle. 


4.2 Contributions to A 

Using the effective vertices from the previous section for calculating the second loop, ignoring sup¬ 
pressed terms proportional to higher powers of rrv^/M 2 (with M a heavy mass) in the numerator and 
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the muon mass in the denominator, we obtain the various contributions to the anomalous magnetic 


M 1 ’ = E 




a m„ 


47r 3 v 2 


• v cQf 


Re 


(vf) R e(yf) + 1 ™{vf) l Hvf) ^ {jtfT 

P i Pi 


moment of the muon. The first two contributions are the well known classical results |6|[7,31-41 

.2 


( 18 ) 


M 2) = E ^§r Mvf) K°»+h- ? {2) ( 


The third contribution simply reads 


K? 


i y fi 


( 20 ) 


( 21 ) 


As for the new contributions, given by the last three sets in Fig. [3j their contributions are given by 


A ( 4 ) _ amlN c \V t b \ 2 r 1 

% “ 32 773 4^^-114) J 0 


dx 


Qtx + Qb (l - x ) 


Re(s d tf) m 2 x(l - x) + R,e(? u c*) m%x(l + x) 


mi 


mi 


M“h± ’ 


-g 


( m t 

m 2 b \ 

Ww 

M 2 V )_ 


( 22 ) 


Aa (5) = __ 

M 64tt3.4u 2 (M4-A4) 


E Re sf 1Zii(TZi2 — iR-is) 


(M|± + A4 - M\) (1 - x) - 4A4 

r i ' 


dx x~ 


/if 2 M 2 0 

Vi 


Af| ± ’ M| ± 


1 , 


M 2 0 

(Y?V 


(23) 


Aa ( 6 ) =__ 

M 64 7T 3 4 (Mfj-t - M\y ) 


y~] Re sf ( 7^*2 - 7 ^: 3 ) \o H+H - / dxx 2 (x - 1 ) 


<5 l 


M 2 0 

_ 

M 2 H± 


-G 


Mg± kk 

M 2 ’ M 2 

iK/ VF iK/ VK 


.(24) 


We can also consider the contribution from a lepton and a neutrino loop by replacing Q t —> 0, m* —> 0, 


Qb —> —1, mfe —>■ mi, Sd —> Q and —> 0 in (22) and where mi is the mass of the considered lepton. 
However, these contributions turn out to be very suppressed due to the smallness of the lepton masses 
and we shall ignore them in our present analysis. The needed loop functions are given by: 




W f 1 dx 2x{l ~ ~ 1 

2 Jo w-s(l-s) 


In 


u 


x{\ — x) 


(25) 


J&Xu) 


- [ dx - 

2 Jo u- x{l 



OJ 


x(l — X) 


(26) 
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(27) 


and 


J^{u) 


1 [ 1 dx hl ( U 

2 J 0 w-x(l-x) \x(l-x) j ’ 




1 f 1 x [3x(4x — 1) + 10]w — x(l — x) ^ ( oj \ 

2 Jo w — x(l — x) l x(l — x) ) ’ 


In 


G(u a : u b ) = 


oj a x + w b (l — x) 
x(l — x) 


x(l — x) — ca a x — w b (l — x) 


(28) 


(29) 


5 Phenomenology 


In the present analysis we neglect possible CP-violating effects; i.e., we consider a CP-conserving 
scalar potential and real alignment parameters The fermionic couplings of the neutral scalar fields 
are then given, in units of the SM Higgs couplings, by 


Uf = cos a + sin a , 


Vd,l ~ 1 Sd,i, 


H 


Dj = — sina + cos a , 


Vu = -is u , 


(30) 


and the couplings to a pair of gauge bosons (6) are simply (k.(/ = g^o w/dhvvi V = W,Z) 


,SM 


Hy = Tin = cosd, At y = 77-21 = —sina, Ky = 77-31 = 0. 


(31) 


We shall separate the phenomenological analysis in two parts. For the first part we will analyse the 
individual contributions from the various A ay factors for different coupling and mass configurations. 
As for the second part we shall sum all these contributions choosing a few relevant scenarios compatible 
with collider and flavour bounds and also with constrains from the oblique parameters. Also, we will 
identify the lightest CP-even Higgs with h and take Mh = 125 GeV for the whole analysis. 


5.1 Individual Aa? contributions 


As we know from global fits to the LHC data, the Yukawa couplings of the discovered scalar boson 
are SM-like, however with quite large experimental errors. The coupling of h to two gauge bosons is 


constrained by | cosd| > 0.8 at 95% CL 43 . Here we shall always take the positive solution, cosd > 0 
(flipping the sign of cos a leads to an equivalent solution with a sign flip of the couplings ?/). Choosing 
the positive solution for cosd, the top Yukawa coupling must also be positive. We shall vary it in the 
range y b G [0.8,1.2]. As we know, at least for now, there is no experimental sensitivity to the relative 
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Figure 6: One-loop scalar contributions to Aas functions their couplings to fermions from h (top-left), H 
(top-right), A (bottom-left) and H ± (bottom-right). 


sign of the down-type or leptonic Yukawas with respect to the up-type Yukawas. Therefore we shall 
be less restrictive with the j/)j t couplings and allow them to vary in the range y^i G [—1-5,1.5]. As 
for the alignment parameters, we will vary them as follows: — 1 < <, u < 1 compatible with all flavour 
constraints and direct charged Higgs searches [43 for a broad range of the charged Higgs mass, and 
—50 < (, c i.l < 50 to safely avoid the non-perturbative regime. We shall also vary yJ in the same 
regions as the parameters (in the limit cosd 1 we obtain y j? = <^). The remaining parameters 
are the couplings of the neutral scalars to a pair of charged Higgses. In order to safely satisfy the 


perturbativity bounds 44 for a broad range of M H ±, we will impose IA^o^+^-1 < 5. 

The one-loop well known contribution from the various scalars are shown in Fig. [6] The contribu¬ 
tion of h is small and positive for the whole considered range of the coupling 1 1 /( \. The contribution 
of H is also positive and, its contribution can be of some significance only for large values of \yj^\ and 
small values of Mu simultaneously. The contribution of the CP-odd scalar is negative and it is only 
relevant for large values of |q | and low values of it mass mass, similar to the previous case. As for the 
charged Higgs contribution, it is always negative and very small, thus irrelevant, at the one loop level. 

The two-loop results are presented next. The contribution of h, associated with a top-quark loop, 
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Figure 7: Contributions to A cffi from h (top-left) and H (top-right) with an associated top-quark loop, and 
from H with an associated bottom-quark (bottom-left) and tau-lepton (bottom-right) loop , as functions of their 
couplings to fermions. 


to A aft is shown in Fig. [7] (top-left). It is positive for y^y^ < 0. The contribution of the same scalar 
h associated with bottom and tau loops is much smaller, of 0(1O~ 13 ) or less for the whole considered 
parameter space, and is not shown here. The contribution of H for different mass configurations 
and for different fermionic loops is also shown in Fig. [7J This contribution is proportional to the 
yj^ coupling which can be large. Thus is turns out to be non-negligible even for the sub-dominant 
bottom-quark and tau-lepton loops. The top-quark loop contribution can be large for all considered 
mass settings as long as yf* is large, and it is positive for yf* y^ < 0, as we can observe in Fig. 0 (top- 
right). The bottom-quark loop contribution can be additionally enhanced by the coupling y^, thus, it 
can overcome the mass suppression. This contribution is positive for yf 1 y^ < 0, see Fig. [?] (bottom- 
left). Similar considerations about the enhancement factor (yf 1 ) 2 can be made for the tau-lepton part, 
however this contribution is always negative, as shown in the bottom-right panel of Fig. [7j 

The contribution of the CP-odd scalar to Aajj 1 is probably the most interesting yet. It has been 


extensively analysed in previous works |6)[7,31 37 . For low values of its mass and large values of ^,1 
it can reach values within or close to the two-sigma region of Aajf p , as it is plotted in Fig. [8j Its value 
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Figure 8: Contributions to A a\}' > from the CP-odd scalar A, associated with a top-quark (left) and bottom-quark 




Figure 9: Contributions to A from gf) = h , H as functions of the product of the couplings y‘f* X v o H + H - for 
various charged Higgs masses (left) and for various configurations of (M#±, Mh) (GeV) (right). 


is positive for q u q < 0 (q q > 0 ) for the top (bottom) quark loop contribution and is always positive 
for the tau loop contribution. This last case is not shown. It is worth mentioning, however that the 
tau loop contribution is somewhat larger than the (absolute value of the) bottom contribution. Even 
if the tau-lepton has a relative mass suppression, the bottom-quark has a charge suppression that is 
in general larger. 

For Aa ^ 2 we only have two possible contributions, from h and H (in the CP-conserving limit the 
vertex AH + H~ vanishes j44j). The contribution of the light scalar h is relatively small for the whole 
considered parameter space, Fig. [ 9 ] (left) and that is due to the fact that y\ l E [—1-5,1.5] whereas yf 
can be much larger. The contribution of H can be quite large depending on the configuration of the 
masses (M H ±, Mh) (GeV). It reaches its largest value for low masses of both M H ± and Mh and large 
values of the product of the couplings y { H A hh+h~ ■ However, even for lower values of the couplings 
but with low masses (or large masses and large couplings) the contribution can be non-negligible. For 
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Figure 10: Contribution to A from H (left) as function of the product of the couplings IZ 21 yf 1 for various 
mass configurations. Two-loop dominant contribution from the top-bottom quark loops to Aaj 1 4 ' 1 (right). 


details see Fig. E (right). 


A (3) 

The next contribution we focus on is A ay,. The contribution from the light scalar is small, of 
(^(lCP 11 ) or less (after subtracting the SM), therefore we can safely neglect it. The H contribution 
however, is non-negligible. It reaches higher values (and it is positive) for low values of Mjj and large 
positive values of the product 7^21 yf 1 (= sin 2 a — q sin a cos a) as it is plotted in Fig. 10 (left). As 


we have already mentioned before, this diagram should not be neglected, as it can introduce sizeable 
effects for some regions of the parameter space. 

Now we move on to the analysis of the charged Higgs contributions of the Barr-Zee type diagrams 
(Fig. §, which is the main goal of this paper. It is obvious from Fig. [lo] (right) that the A 
contribution is non-negligible for a large region of the parameter space, except for very small values 
of the product For a charged Higgs with a low mass, say 90 GeV, and large negative values of 

q q u this contribution alone can explain around 35 % of the measured discrepancy. This looks very 
appealing, because with the exception of a very light CP-odd scalar, the previous contributions cannot 


reach such large values. For the plot shown in Fig. 10 (right) we have chosen = 0. However, a 
variation of Sd in its allowed interval [—50, 50] only produces a shift in the plotted values of order 10 -12 
or less. This is obviously due to a relative suppression factor m 2 /m 2 and therefore this contribution 
can be safely ignored. 

Last, contributions A affl and Aa[f' 1 are shown in Fig. |ll|and Fig. 12 They are a little bit smaller, 
however they can reach values up to 10~ 10 . Again this happens, for small mass configurations and 


large values of the corresponding couplings. We can see in Fig. 11 that both h and H contributions can 
be very similar, however, they cannot be simultaneously positive (if the product of the three couplings 
q TZj 1 TZ ,;2 is chosen positive for one scalar, for the other must necessarily be negative). On the other 
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SV ^11 ^12 Q ^21 ^22 

Figure 11: Contributions to A from tp® = h,H as functions of the product of the couplings q IZu Ha for 
various charged Higgs masses (left) and for various configurations of (Mh, M h ±) (GeV) (right). 




Figure 12: Contributions to A from pf) = h, H as functions of the product of the couplings <Tz Ha ^ v o H + H - 
for various charged Higgs masses (left) and for various configurations of (Mh, M h ±) (GeV) (right). 

hand, both h and H contributions from A affi can be simultaneously positive, and of similar value. 
Thus, when summed up they can play an important role in the total value of A a^. 

We have proven thus, that these new Barr-Zee contributions must not be ignored, as they might 
sizeably modify the theoretical prediction for this observable within the 2HDM framework. 


5.2 Total contribution to (g — 2) M 

Thus, we have seen that the dominant contributions of the new Bar-Zee type diagrams come from 
the mechanisms (3) (Fig. [2]) and (4) (Fig. [3]). All the other new contributions are sub-dominant. 
Now, it is interesting to put all these results together, and show the total effect on Afor a few 


relevant scenarios. In Fig. 13 (left panel) we show Aas a function of q for positive values of 
this coupling and for a few scenarios given by cosd = 0.9, s u = —0.8, q = —20, Mg = 125 GeV, 
A h.H+H- = 0 ; A hH+H - = — 5. The masses (in GeV) of the remaining scalars are chosen the following 
way: Mh = M H ± = Ma = 250 (lower orange curve), 150 (middle blue curve), Mh = M H ± = 150 and 
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Figure 13: Total Aa M contribution as a function of q /or different coupling and mass configurations. 


Ma = 50 (upper green curve). Similar to the previous case, but this time for negative values of q, in 
the right panel we have chosen the following parameter configuration: cos a = 0.9, s u = 0.8, Sd = 2, 
Mh = 125 GeV, A hH+H- = 0; A hH+H- = 5 and Mu = M#± = 250 GeV and Ma = 40 GeV (upper 
green curve) or Mu = M H ± = 350 GeV and Ma = 50 GeV (lower orange curve). As expected, from 
the analysis of the various A af individual contributions, one obtains a significant contribution for 
low masses of the scalars (especially for low Ma) and large couplings. We can also observe that in 
some cases we do not need the maximum allowed value of |q| in order to reach the two-sigma region 
of Aa value around |q| ~ 30 might just be enough. 


6 Conclusions 


It is a common belief that only a restrained number of diagrams, namely (1) and (2) from Fig. [2] 
can significantly contribute to Ain 2HDMs and in most of the previous analyses @01 3lJ-37j, a 
CP-odd scalar in the low-mass range is enough to explain, or reduce, the discrepancy between theory 
and experiment. In this work we have shown that the extra degrees of freedom of the A2HDM given 
by the Sf parameters, can also explain this discrepancy in some region of the parameter space, and if 
not, they can significantly reduce it in most cases. We have also seen that the W loop contribution 
associated with a heavy scalar H (diagram (3) from Fig. [ 2 ]) can bring important contributions even 
if it has a global suppression factor 72.21- This contribution is positive for negative values of q. The 
most interesting case is, however, the fermionic loop contribution (diagrams (4) from Fig. [ 3 ]) with the 
dominant part given by the top-quark. The last two diagrams (5) and (6) are also interesting, as 
they can sum up to an 0(10%) of the total contribution. Also, we have seen that not all of these 
new contributions can be made simultaneously positive, however the total Ais positive for most 
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parameter configurations. 

A highly interesting scenario, that we defer for future work, is to consider CP-violating effects. 
The imaginary part of the parameters of the potential and especially of the Yukawa sector might be 
able to bring somewhat sizeable effects. 


A WW 7 effective vertex contribution to (g — 2) M 


In this section we present the explicit calculation of the contributions from Fig. 0 (B) to (g - 2) m . 
The 2HDM contributions to the one-loop WW'y effective vertex are shown in Fig. |TT| where last 
diagram stands for the one-loop renormalization counter-term. For this calculation we have followed 
the renormalization prescription described in [61 . Following this prescription one does not need to 


renormalize the gauge-fixing Lagrangian. Thus, we simply worked in the Feynman gauge 44 . Working 


in this gauge, one also needs to take into account IFG^y (Fig. 15) and G^G^y effective vertices. The 
last set (G^G^y) will give rise to contributions to the anomalous magnetic moment that will have a 
relative suppression factor of (just as in case (A) of Fig. [ 4 ] for the 7 effective vertex), 

and therefore will not be taken into account. 

The one-loop counterterms for the needed IF IF 7 and IFG 7 vertices are given by 


* r T = i e S w , iT? = ie g ^ - (6 W + 6 G ± + S M ), (32) 

where ieT ppv is the tree-level WW'y vertex and where we have defined the G f± , W p and renor- 



Figure 14: One-loop contributions to the WW'y effective vertex. The last diagram stands for the one-loop 
counter-term. 
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Figure 15: One-loop contributions to the W 7 effective vertex. The last diagram stands for the one-loop 
counter-term. 


malization constants as 


Z\v = 1 + 5w , 


Z G ± — 1 + 5 g ± , 


Zm — 1 + 5m ■ 


(33) 


The needed W and G± self-energy diagrams needed for the calculation of these counter-terms are 


shown in Fig. 16 As we can see, no tadpole diagrams are present. At one-loop level, using the renor¬ 


malization prescription from 61 , tadpole diagrams do not contribute to the W mass renormalization. 


On the other hand, they do not contribute to the wave-function renormalization either as they do not 
generate any four-momentum dependence. Thus, for our present calculation we need not to worry 
about tadpoles. 

One last technical issue is the W — mixing that occurs at one-loop level. The gauge fixing 
Lagrangian cancels exactly the tree-level mixing between the gauge and Goldstone bosons generated 
by the covariant derivatives. This mixed term, when renormalizing the Lagrangian is in fact, counter¬ 
term for the W — G f± self-energies, as it is nicely explained in 1 6 l]. For this calculation, however, 
we don’t need to worry about this mixture. As we are going to ignore the propagator corrections, 
and these corrections are related to the W — Cr mixing through the Ward identities (for example 
the doubly contacted identity shown diagrammatically in Fig. we are also going to ignore the 
one-loop mixing in order to preserve these identities. 

Using the on-shell scheme, working in D = 4 + 2e dimensions (e < 0), the expression for 5w reads 

5w = 4 / + 5 w + 44 with: 


= IK<2 + n a \ 


2 M 


2 e 


(47r) 


h + 


dx 2 x'(l — x ) In ■ 




w> 




(34) 
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Figure 16: One-loop W self energy diagrams needed for the vector boson wave-function and mass renormaliza¬ 
tion, and G ± self-energy diagrams needed for the charged Goldstone boson wave-function renormalization. 
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Figure 17: One-loop doubly contracted Ward identity. The grey circles stand for the one-loop self energies. 


A2) _ M'w 
°W - v 2 


x(3) _ 4 

°W - v 2 


tj^(h + L 41 21(1 “ 1)111 

? K?i w? L 


T 2 


" 1 , x(l — x) 

dx ^m 


(35) 


(36) 


in agreement with [621. The wave function renormalization counter-term for the charged Goldstone 
boson is given by S G ± = with: 


/ 0 a 2 ( M w) ’ 


c(i) _ 1 \V 12 Mr± ^ [ l , x(l-x) 

^G± (47t) 2 1^*2 + ^. 3 1 - / dx 

°G± ( 47r )2 


M 4 n ,1 


__ 1V1 0 /* 


dx 


x(l — x) 


0 a 2 ( M w) ’ 


(37) 

(38) 


,2e 


4± = “7T7W [ f + l + l dx ( 3x2 - 6x + 4 ) ln ° ^2 W) 


( 47r) 2 i ; 2 u L e ' 6 

v i 

2 ^2 


+ 


dX a 2 (M 2 ) ( M wfi x2 “ 8x + 6 ) + 2xM ^a) 


(39) 


Last, the VF mass counter-term is given by 5 m = <5^ + <5^ + 4” with: 
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2 e 




V 2 (47r) 2 


E 1^*2 + ^"t3| 2 ( + - l) { M H± + M*o - ^M 2 V ) 




^( 2 ) _ 1 ^ 

— 9 

V z 


+ / cfa2o 2 (Ma,)ln^^ 

do I- 1 


r-l 


° 2 (Myy ) 


2 E -^1 ( g “ *) «o + - Mw) + jf dx 2 a 2 (Mw) In — 


r(3) 4M j+ r>2 

M v 2 (4vr) 2 ^ ?:l 


1 f 1 

t + / dx In 


J (M 2 


W/ 


/T 


( 4 ) _ 2M 2 h± y 2e 

M v 2 (4vr) 2 


M 2 


1 + ln^fli-1 

e /i z 


x(5) _ 

- (47r)2 


E 


Mjrl 

x 2 I. e 


M 2 o 

+ In —- i 
P 


(40) 

(41) 

(42) 

(43) 

(44) 


Here we have defined 1/e = 1/e + 7 e — ln( 47 r). The functions a 2 (p 2 ) and a 2 (p 2 ) are given by: 


<*V) 

= —p 2 x(l — x) 

+ M 2 o X 

+ Mft± (1 — x) , 

(45) 

d 2 {p 2 ) 

= — p 2 x(l — x) 

+ M 2 o x 

+ M\v (1 ~~ x) . 

(46) 


Now, we move on and present the expressions for the one-loop WW'y effective vertices from Fig. [14 
The considered kinematics and the assigned Lorentz indices for this process are W + (k—q, p)+ r y(q, p) —5 
W + {k,v). Discarding all terms proportional to q p , the first and second diagrams give: 


i TZ U = i 


e M 2 v 2 £ 


(!) ‘(4t r) 2 v 2 


p 2e E l 7 ^ 2 + K l5 \ 2 \- + f 1 dx [ l dy 2(1 - x) In ^ 

L rfe Jo Jo P 


+ f 1 dx C dy 
Jo Jo 

-1 /•! 


jpfU' 

J (b) 


k 2 — M 2 — 2 y k ■ q . 


(47) 


> rgr = E R *i [ - M r "“' + / d \[/y •'w 2 <1 - *>I 


+ f 1 dx f 1 dy 
Jo Jo 


T ppu 

J (b) 


k 2 — M 2 — 2 y k ■ q - 

Again, r p/iI/ is the tree-level vertex function and it is given by 

T ppv = <?^(-fc - g) p + g pp { 2 q - k) v + g vp ( 2 k - q) p . 

The sum of diagrams (3), (4) and (5) gives 

E R ? i / 

do 


• pP/ 2 i 2 _ ■ 6 

(3+4+5) (4vr ) 2 X 2 


1 r 1 1 -2g pv Q p -2g pp J v + AJ P ^ 

dx / dw---=—- 

do * k 2 — M 2 -2yk-q 


(48) 


(49) 


(50) 
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With diagram (6) we have to be specially careful. Its explicit expression reads 
e Mb 


• v p^ = • 

( 6 ) (4vr) 2 v 2 


^ - <f) g 2e Y, 1^2 + ^3| 2 dx (2x - 1) (t + In . (51) 


Integrating over x, the pole and the /i-dependence vanish. We are left with a logarithm that depends 
on the four momentum and that we need to integrate in the second loop. Using the expansion (d <C 1) 

1 


lnA= |(a*-i) + 0 ( 6 ) , 


(52) 


we can write the previous expression as 
e Ml 


j rP^ = j 

( 6 ) ( 47 t ) 2 V 


f m.+<*«r jf * ■ (53) 


and use the Feynman parametrization 

1 

a a tv r\ 

1 0 

X 


p/Q _ X) r 1 rl-xi rl-x\—X2 

= TMf / ^ / 

1 


(54) 


(xiAi + x 2 A 2 + x 3 A 3 + (1 - Xl - x 2 - x 3 )yf 4 ) 3 5 
in order to solve the second loop (taking the limit <5 —>• 0 at the end of the calculation). We obtain a 
similar expression for diagram (7): 

e M 2 


i = i- 


_ _ I w a ^ (k P_ a P ) t}) S . V77 2 f 1 dx ( 2x ~ l ) x5 ( l ~ x ) 5 (55) 

( 7 ) ‘(4vr) 2 v 2 9 Q) 5 ^ Uil J 0 d [k 2 - M 2 -2k- q)~ 5 ' (55) 

Contributions (8) and (9) vanish as their expressions are terms proportional to q p . Finally, diagrams 
(10) and (11) read 

^ ^|fti 2 + ^ 3 | 2 J dx (2x - 1) (1 + ln(|) » (56) 

s'Vc-i;*., I'dxi 2*-1)(I+!..)§) , (57) 


which can be treated exactly as diagrams (6) and (7). The previously introduced tensorial functions 
are given by: 

J p $ = 1 - 2 x)k v + 2 y(x - 1 )q u ) + g^(( 1 - 2 x)k p + (2(x - l)y + l)<f) - 2x ^ ^ , (58) 

= -2P((2x - l)*" - 2y(x - l)<f) ((1 - 2x)F + (2(x - 1 )y + 1 )q p ) , (59) 

= g pp ((xy — y + 2)(/ t/ — xAT) — g pu (xk p + q p (y — xy + 1)) + 2x /c p , (60) 

and, 


Q p = (1 - 2x) + q p (2xy - 2y + 1) 


J 13 = AT (1 - 2x) + q v 2 y{x - 1) . 


(61) 


22 



The scalar functions are given by: 


«x = -®(! - x ){k 2 - Ml - 2 y k ■ q) , 
bl = — x(l — x)(k 2 — M 2 — 2 k ■ q) , 
cl = -x(l -x)(fc 2 - Ml) , 


a l = ~ x i l ~ x )( k2 ~ Ml -2 yk-q), 

bl = —x{l — x)(k 2 — Ml — 2 k ■ q) , (62) 

c 2 x = - x (\-x)(k 2 - Ml) , 


with 


, r9 K° m 2 h± 

Ml = —+ -ML_ 
1 — x x 


2 

2 _ A , 


Mi = 


M l* Ml 


1 — x 


+ 


(63) 


Next we present the G^W'y effective vertices from Fig. 15 The kinematics and Lorentz indices are 
given by G + {k — q) + 7 (q, n) —> W + (k, v). Thus, the one-loop expressions are: 


i rg = — ^ 4^)2 Mw H 2e ^2 |^»2 + 777* 3 | 
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The tensorial functions are given by: 


\ + ^ dx In 4 ) • 
e Jo ^ y 


(71) 


= k» ((2x - 1) F - 2y(x - 1) q v ) , (72) 

G MI/ = g^ u [k 2 x(x — 2) — 2(x — 1 )(xy — y — 1 )k ■ q ) + /A((x — l)(xy + 2y — 4 )q u — x(x — 2)k v ). (73) 


All other functions are the same as previously. Note, that for the previous expressions of the one-loop 
effective vertices, we have maintained the k■ q structure in the denominator (in contrast to the H±Wj 
effective vertices) because here, in some cases, this structure does contribute to the final result. 

Inserting all these expressions into the second loop we finally obtain the expression for the total 
contribution to the anomalous magnetic moment of the muon. Subtracting the SM contributions we 
have 


A " 


dx f 'y ] A — Asm + ^ ' \T^i 2 + *7431 2 B + C 


128 tH 4 v 2 J 0 
with the functions A, B and C given by: 

„ 7 a 2 (M 2 ) (2x 2 - 3x + 2) A/ 2 0 6 (x - 1) M 2 + (-12x 2 + 30x - 55) M 2 

^-sjsrrjg 1 +-' 6 (a4,-a5)-- + 
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All the functions that carry a SM sub-index are obtained from the original ones by replacing M o with 

f i 

everywhere, where is the mass of the SM Higgs. The numerical values that we obtain for this 
contribution (for M H ^ H ± < 500 GeV) are typically of 0( 10~ n ) both positive or negative, which is 
two orders of magnitude below Aa)) xp , therefore we shall not take it into account in this analysis. 
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